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b-H ■ Abstract. It is shown how to define difference equations on particular lattices {x n }, 

\ n £ Z, made of values of an elliptic function at a sequence of arguments in arithmetic 

■ progression (elliptic lattice). Solutions to special difference equations have remarkable simple 

CNJ \ intcrpolatory expansions. Only linear difference equations of first order are considered here. 
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1 Difference equations on elliptic lattices 

m : 

O ■ 1.1 The difference operator 

OO ■ 

We consider functional equations involving the difference operator 

o; (vf)(x)- m 

, Most instances [26] are (f(x),ifj(x)) = (x,x + h), or the more symmetric (a; — h/2,x + h/2), 

or also (x,qx) in q-difference equations [131 CE1 E! • Recently, more complicated forms (r(x) — 
^Jsjx), r(x) + ^/sjx)) have appeared PQEl ESI IIZl E21 EZ31 C2ZI EBl > where r and s are rational 
5^ , functions. 

This latter trend will be examined here: we need, for each x, two values f(ip(x)) and f(ip(x)) 
for /. A first-order difference equation is 

T (x, f(<p(x)), f(ip(x))) = 0, or f(<p(x)) - /(>(*)) = Q (x, f(<p(x)), 

if we want to emphasize the difference of /. There is of course some freedom in this latter 
writing. Only symmetric forms in (p and ip will be considered here: 

(Vf)(x) = ^(x,f(tp(x)),f^(x))), 

where T> is the divided difference operator ([1]) and where & is a symmetric function of its two 
last arguments. 



*This paper is a contribution to the Proceedings of the Workshop "Elliptic Integrable Systems, Isomonodromy 
Problems, and Hypergeometric Functions" (July 21-25, 2008, MPIM, Bonn, Germany). The full collection is 
available at |http://www. emis.de/journals/SIGMA/Elliptic- Integrable-Systems.html] 
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For instance, a linear difference equation of first order may be written as 
a(x)f(<p(x)) + b{x)f{i>{x)) + c(x) = 0, 
as well as 

a(x)(Vf)(x) = /3(x)\f(<p(x)) + f(i/>(x))] + 7 (x), 

with a(x) = [b(x) — a(x)][ip(x) — tp(x)]/2, 0(x) = — [a(x) + b(x)]/2, and 7(2;) = —c(x). 

The simplest choice for ip and ip is to take the two determinations of an algebraic function of 
degree 2, i.e., the two y-roots of 

F(x,y) = X {x) + X^y + X 2 (x)y 2 = 0, (2a) 

where Xq, X\, and X2 are rational functions. 

Note that the sum and the product of (p and ip are the rational functions 

tp + ^ = -X 1 /X 2 , ipip = X /X 2 . (2b) 



1.2 The corresponding lattice, or grid 

Difference equations must allow the recovery of / on a whole set of points. An initial-value 
problem for a first order difference equation starts with a value for f(yo) at x = xo, where yo is 
one root of ([2a|) at x = xq. The difference equation at x = xq relates then f(yo) to f(yi), where 
yi is the second root of (|2"a|) at xq. We need x± such that y± is one of the two roots of (|2"a|) at x±, 
so for one of the roots of F(x, y±) = which is not xq. Here again, the simplest case is when F 
is of degree 2 in x: 

F(x,y) = Y (y) +Y 1 (y)x + Y 2 (y)x 2 = 0. (2c) 

Both forms (|2ap and (j2c|) hold simultaneously when F is biquadratic: 

2 2 

i=0 j=0 

The construction where successive points on the curve F(x,y) = are (x n ,y n ), (x n ,y n+ i), 
(x n+ i,y n+ i), is called "T-algorithm" in [34\ Theorem 6], see also the Fritz John's algorithm 
in [UE1E]. The sequence {x n } is then an instance of elliptic lattice, or grid. 

Of course, the sequence {y n } is elliptic too, x n and y n have elliptic functions representations 

x n = £1 (t + nh), y n = £ 2 (t + nh) , (4) 

where (x = £i(t),y = £2^)) is a parametric representation of the biquadratic curve F(x,y) = 
with the F of ©. 

Note that the names of the x- and y-lattices are sometimes inverted, as in |344 equation (1.2)] 
As y n and y n +i are the two roots in t of F(x n ,t) = Xq(x u ) + Xi{x n )t + X 2 (x n )t 2 = 0, useful 
identities are 

X\{x n ) _ X ( 

X 2 {x n ) X 2 (x n ) 

from (|2b[) . and the direct formula 

, -X 1 (x n )± v/pg) 

y n and y n+ i = v , 

2X 2 {x n ) 
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where 

P = X\- 4X X 2 

is a polynomial of degree 4. 

Also, as x n+ i and x n are the two roots in t of F(t,y n+ i) = 0, 

_ Yi{yn+i) _ Y (y n+1 ) 

X n + Xn+i — . . , X n X n +i — . . . 

Y 2 {yn+i) Y 2 {y n+1 ) 

As the operators considered here are symmetric in <p(x) and ip(x), we do not need to define 
precisely what tp and ip are, i.e., we only need to know the pair (ip, ip), and not the ordered pair. 
However, once a starting point (xo,yo) is chosen, it will be convenient to define p(x n ) = y n and 
ip(xn) = Vn+i, neZ. 

Special cases. We already encountered the usual difference operators (p(x),ip(x)) = (x, x+h) 
or [x — h,x) or (x — h/2,x + h/2) corresponding to X 2 (x) = 1, X\ of degree 1, Xq of degree 2 
with P = Xf — 4XqX 2 of degree 0. For the geometric difference operator, P is the square of 
a first degree polynomial. For the Askey- Wilson operator [TJ [21 [T21 121 E21 E3] , P is an arbitrary 
second degree polynomial. 

The formulas for the sequences x n and y n are in these three special cases 

(x n ,y n ) = (x + nh,yo+nh); (a + bq n ,u + vq n ); 
(a + bq n + cq~ n ,u + vq n + wq~ n ). 



1.3 Difference of a rational function 



From (|2bj) . when the divided difference operator D of ([T]) is applied to a rational function, the 
result is still a rational function. 

The difference operator applied to a simple rational function is of special interest. 

Let f(x) = -^j, then 



v 1 



x — A ip(x) — tp(x) 



1 1 



ip(x) - A ip(x) - A 
X 2 (x) 



Wx) - A)(p(x) - A) 



X (x) + AX 1 (x) + A 2 X 2 {xY 

and let {{x' n ,y' n ), ( x' n , y'n+i)} be the elliptic sequence on the biquadratic curve F(x,y) — such 
that y' Q = A, then 

v-L- = Mx) (5) 

x-A Y 2 {A){x - x' ){x - x'^Y K ' 

as the denominator is F(x,A), and the two x-roots of F(x,A) = F(x,y' Q ) = are x' and x'_±, 
from the opening discussion of Section 11.21 

The T> operator applied to a general rational function yields a rational function with the 
factor X 2 . It seems sometimes fit to define a difference operator as our T> divided by X 2 , as by 
V.P. Spiridonov and A.S. Zhedanov in Section 6 of [32]. See also Section 2 of [34] . 

A general rational function is generically a sum of simple rational functions of type ©, say, 
l/(x — A), l/(x — B), etc. The difference has poles at x' and x'_ 1 , also at x'q and x"_ x if 
B = t/q, etc., so that the degree of T>f is usually twice the degree of /. However, the difference 
of a rational function of denominator (x — y' )(x — y[) ■ ■ ■ (x — y' n ), T>f has no other poles than 
x'-i, x' , . . . , x' n . This is also discussed in 



4 



A. P. Magnus 



So, let {(x n ,y n ), (x n ,y n+ i)} be a first elliptic sequence on the biquadratic curve F(x,y) = 0, 
and {(x' n ,y' n ), (x' n ,y' n+l )} be another elliptic sequence on the same curve. The two sequences 
have the same formula (J3J), but with different starting values to an d t' . 

Now, let 

v / n (x-x )---(x- x n -x) (x-y )---(x- y n -\) 

x n(x) = — r j— and y n (x) = — r — . 

(x-x x )---{x-x' n ) (x-y x )---(x-y' n ) 

See that 

X n -\{x) 



vy n (x) = c n x 2 ( x ) 



(x -x' ){x -x' n y 



Indeed, (ip(x)-y )(ip(x)-yi) ■ ■ ■ (tp(x)-y n -i) and (ip(x)-y )(ip(x)-yi) ■ ■ ■ (ip(x)-y n -i) both 
vanish at x = xq, x\, . . . , x n - 2 ; (y>(x) — y[)(ip(x) — y' 2 ) ■ ■ ■ {fix) — y' n ) vanishes at x = x[, . . . , x' n , 
whereas (i/j(x) — y' 1 )(ip(x) — y' 2 ) ■ ■ ■ (ip(x) — y' n ) vanishes at x = x' , . . . , x' n _ x . 

Simple fractions give 

v _ 1 X 2 (x) 



x-y'k Y ^v' k )( x - x ' k -i)( x - x ' k ) : 

as seen earlier in (|5|). 

The constant C n is found through particular values of x, either x-i, where y n {ip{x)) = but 
y n ((p(x)) 7^ 0, or x n -i, where y n {ip(x)) = but y n (ip(x)) / 0: 

c = y n (<p(x-i) = y-i)(a_i - x )(x_i - x' n ) 
(j/o - y-i)X 2 {x-i)X n - 1 {x- 1 ) 

c = 3^W^n-l) = y n )(x n -l - x' Q )(x n -i - x' n ) 

(y n - y n -i)X 2 ( Xn— \)X n — \{x n — x) 
(of course, Co = 0). Or through residues at x' , where y n (ijj(x)) = oo, or x' n where y n ((p(x)) = oo, 

c (V'l -VP)--- (V'l - Vn-l) x 0- x n , fi v 

n ^(^-^•••(yi-<)^-^(x )^- 1 K)' 1 j 

c = (y'n -VP)--- (y'n ~ Vn-l) x' n - x ' Q 

n ~ (^-^•••(^-^ 1 )^(^l-^ 2 (<)^l(<)- 

We shall also need the operator M. defined as 

(Mf)(x) = [f(<p(x)) + f^(x))]/2, 

which sends rational functions to rational functions too, usually of double degree, but without 
particular factor. 

With this operator M, 

2(My )(x) - ~ vMpW - yi) ■ ■ ■ (v(x) - vn-i) 



(f(x) - y[)(f(x) -y' 2 )--- (<p(x) - y' n ) 
(ip(x) - y )(ip(x) - yi) ■ ■ ■ (tp(x) - y n -i) 



+ 



(ip(x) - y[)(ip(x) -y 2 )--- (ip(x) - y' n ) 



2D n (x) {x - X0 ][ X - Xl ]''' {x - Xn ;" ) = 2D n (x)- Xn ~ l{x) 



(x - x' )(x - x[) ■ ■ ■ (x - x' n ) n {x - x' ){x - x' n Y 

where D n is a polynomial of degree 2. 
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Interesting values are found at the same point as in ([6]): 

i?n(x-i) = - W -f ~ y - l} , (7a) 
A^-i) = c ^( Xn ^yn-y n -,) ^ (?b) 
DM) = CnX 2 ( X ' M-y' ) ^ (7c) 

DM = y M<) ^-^\ (7d) 
when ?7, > 0. Of course, -Do = 1- 

2 Elliptic hypergeometric expansions 

Let us consider expansions of the form 

oo 

fe=0 j 

where z£ is a combination ajX^ +bj or a,jy£ +bj, {. . . (x^\y^), (x^\ yj^), • • • } being elliptic 
lattices, or grids, related to a biquadratic curve ([3]), the same curve for each j. 

We certainly recover at least a special case of current elliptic hypergeometric expansions, as 
introduced in [U [5j EH EH US . 

2.1 Rational interpolatory elliptic expansions 

Rational interpolants of some function / at yo, y±, . . . , with poles at j/j,^, • • • , are successive 
sums 

,/ s . a-yo x-y (x - y )(x - yi) 

co = j{yo), c + ci T , c + ci r + c 27 7T7 7\> 

x-yi x-y[ {x - y[){x - y' 2 ) 

oo 

5> fc 3> fc (x). ( 8 ) 

k=0 

If, by chance, c& shows a similar form of ratio of products, we see special cases of hyper- 
geometric expansions! This will happen when one expands solutions of difference equations 
which are simple enough. Putting the expansion in the difference equation results in recurrence 
relations for cjt, and we look for cases when this recurrence relation only involves two terms c& 
and Cfc+i. 

2.2 Linear 1 st order difference equations 

a(x)(Vf)(x) = c(x)(Mf)(x) + d(x) (9) 

Where is bl The full flexibility of first order difference equations is achieved with the Riccati 
form 



a(x)(Vf)(x) = b(x)f(ip(x))f(^x)) + c(x)[f(ip(x)) + f^(x))] + d{x) 

but only linear equations will be considered here. However, Q already allows elliptic exponen- 
tials (c(x) = a(x)) or logarithms (c(x) = 0). 
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We now try to expand a solution to ([9]) as an interpolatory series. If the initial condition 
is f(yo) at x = xq, the difference equation allows to find 

j,s = [a(xo)/ (Vi ~ Vo) + c(x )/2]f(y ) + d(x ) 
a(xo)/(yx - yo) - c(xq)/2 

This works fine if no division by zero is encountered. Let us call x' Q one of the roots of the 
algebraic equation 

a(x) c(x) . , , 

S-L = o, at x = x' (10) 



ip(x) — ip(x 



and let, as usual, ip(x' ) = y[, tp(x' Q ) = y' . This shows that y[ is a singular point of /, as trying 
to compute f(y[) from f(y' ) requires a division by zero. Then y' 2 ,y' 3 , ... are poles as well. 
That's why the expansion in ([5]) starts with poles at y^, y' 2 , ■ ■ ■ We also see that such expansions 
represent meromorphic functions with a natural boundary made of poles. At least, if the poles 
are spread on a curve, this will be discussed in Section El 

We also manage to have the initial value f(yo) completely determined by the equation, i.e., 
independent of so, considering 

a(x-i)/ (y - y-x) - c(x-i)/2 

we ask rr_i to be a root of 

a(x) c(x) , . 

+ = 0, at x = x-i. (11) 



i/j(x) — <p(x) 2 

Finally, we shall need the polynomials c and d to be of degree 3, with X 2 as factor: 

c(x) = ((3x + i)X 2 (x), d(x) = (5x + e)X 2 (x). (12) 

We now have enough information for understanding the 

Theorem 1. The difference equation Q on the elliptic lattice F(x n ,y n ) = of (|2a[) (f3|) . where 
a, c, and d are polynomials of degree ^ 3, X 2 being a factor of c and d as in ()12p . has a solution 
with the formal expansion ([8]), where x~\ is a root of (II ip and x' is a root of (|10p . with 

j,, \ d(x-i) d(x-i) + e 

co = f(yo - 



Cl 



a(x-i)/(y -y-i) - c(x_i)/2 c(rc_i) /3x_i+7 
(5 + pc )(x - x[) ( 7 <5 - 0e)( yi - y[)X 2 {x' ) 



Ci(a(x ) - c(x )(yi - yo)/2) (yi - y' )(x - x' )[a(x ) - c(x )(yi - yo)/2] ' 
and when n~^\, 

= Ci x' n -x n ^ r ^ a(x' k ) + c(x' k )(y' k+1 -y' k )/2(x k -x_ 1 )(x k -x' ) 
x[-x C n a (x k ) - c(x k )(y k+1 - y k )/2 (x' k - x-i)(x' k - x' ) 

Ci , _ (y-i -y[)--- (y-x - y' n _ l )X 2 (x^i){x^ 1 - x ) ■ ■ ■ (x_i - x n _ 2 ) 
X x' x -xq n n ~ l (y-i- y%) ■ ■ ■ (y-x- y n - 2 )(x-t- x' ) ■ ■ ■ (x-x- x' n ) 

x yr Q K) + c(x' k )(y' k+1 -y' k )/2 ( Xk - x _ 1 )( Xk - x > ) 
k } a(x k ) - c(x k )(y k+1 - y k )/2 (x' k - x-x)(x' k - x' ) ' 
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Proof. Put the expansion (JSj) in 

oo 

d(x) = a{x)Vf{x) - c(x)Mf{x) = c n [a(x)Vy n (x) - c(x){My n {x)} 



00 X i(x) 

= -c c(x) + ^ c n [a{x)C n X 2 {x) - c(x)D n (x)\ 



(x - x' Q ){x - x' n )' 



The polynomial a(x)C n X 2 (x) — c(x)D n (x) = [a(x)C n — (fix + j) D n (x)]X 2 (x) already has X2 as 
factor from <\12\) . A factor of degree ^ 3 remains. Complete factoring follows: 
at X-i, from ((Taj) and (jTTj) . 



a(x)C n X 2 (x) - c(x)D n (x) = C n X 2 (x^i)[a(x^ 1 ) + (y - j/_i)c(z_i)/2] = 0; 



at x' , from (ITcj) and (fTUj) . 

a(x)C n X 2 (x) - c(x)D n {x) = C n X 2 (x' )[a(x' ) - (y[ - y' )c(x' Q )/2] = 0. 
Therefore we have three factors of first degree 

a{x)C n X 2 (x) - c(x)D n (x) = X 2 (x)(x - - x' )[£ n (x - x n -i) + i]n(x - x' n )], 

where from (I7dl) 



, a{x' n )C n X 2 {x' n ) - c{x' n )D n {x' n ) a{x' n ) + c{x' n ){y' n+l - y' n )/2 



X 2 (x' n )(x' n x—i)(x' n x'o)(x' n x n ^i) (x' n X—l){x' n Xq)(x^ X n —\) 

and from (|7bj) 

a(x n -i)C n X 2 (x n _i) - c(x n _i)D 
X 2 (x n _i)(x n „i - x_i)(x n _i - x' )(x n ^i - x' n ) 
_ c a(xn-i) - c(x n -i)(y n - y n -i)/2 

( x n—l ~ x — l)( x n— 1 — x o){ x n—l ~ x n) 

Next, 

= a{x)Vf(x) - c(x)Mf(x) - d(x) 

00 X (x) 

= -c c(x) - d(x) + V] c n [a{x)C n X 2 (x) - c{x)D n (x)\ ^ — 

1 [x - x Q )[x - af n ) 

= — coc(rr) — d(x) 

+ ]T c n X 2 (x) [£ n (x - Xn -i) + rj n {x - x'J] ^ ^-tfj ^""^ 

(x - x_i)(x - x ) ■ ■ ■ (x - x n - 2 )(x - X n ^i) 
(x - x[) ■ ■ ■ (x - x' n ) 



-c c(x) - d(x) + X 2 (x) 2J c n i 



^ ( x _ x_i)(x - Xq) ■ ■ ■ (x - x n _ 2 ) 

-c c(x) - (f(x) + ciX 2 (x)?7i(x - x_i) 

(x - x_i)(x - x ) • • • (x - x n _ 2 )(x - x n _i) 



{x - x[) ■ ■ ■ {x - x' n ) 
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(x - X-i)X 2 {x) 



-c P -5 + c\7)\ + y^(c + C n+ ir] n+ i)X n ( 



X 2 is a factor everywhere, from (|12p . so 

= -co(flx + 7) - (Sx + e) + C1C1 



00 

+ y^(Cn^n + Cn + ir] n+1 )X n (x), 
1 



a(x ) - c(x )(yi - yo)/2 



Xq — X 



(x - X_l) 



<i(x_ 



1, 



d(x_ 



as 



ci 



Cn+l 
Cn 



a(x-i)/(yo ~ y-x) - c(x_i)/2 
(5 + /3c )(x - x[) 



(5x-i + e)/(/3x_i + 7), 



c(x„i) 
(7<5-/3e)(yi-yi)X 2 (x ) 



Ci(a(x ) - c(x )(yi - yo)/2) (yi - y' )(x - x' )(o(x ) - c(x )(yi - yo)/2) ' 

tin _ C n a{x' n ) + c(x' n )(y' n+1 - y^)/2 (x w - x_i)(x n - x )(x n - x' n+1 ) 
Cn+i a(x n ) - c(x n )(y n+ i - y n )/2 (x' n - x_i)(x' n - x' )(x' n - x ra -i) : 
' l a(x' k ) + c(x' k )(y' k+l -y' k )/2 



Vn+X 



X n -X) 



n 



-A'„(x_i)A'„(x ). 



C n LL a(xfc) + c(x fc )(y fc+ i - y k )/2 

The formula (|13p achieves a construction of hypergeometric type, as each term is a product 
of values of elliptic functions with arguments in arithmetic progression. The exact order of each 
term, i.e., the number of zeros and poles in a minimal parallelogram, is not obvious [33]. Of 
course, a factor like, say, x_i — x k is an elliptic function of order 2 of to + kh from (|4|). The 
same order holds for the ratio 

x_i — x k x_! — £i(t + kh) 

y-x-yk y-x - £2(^0 + kh) ' 

as zeros of the numerator and the denominator cancel each other. 

Similar effects probably hold in other ratios encountered in (|13p . such as 

a(x fc ) - c(x k )(y k+1 - y k )/2 
(x fc - x_i)(x fc - x' ) 

but it is not clear if more can be obtained by keeping elementary means, or if more elliptic 
function machinery (theta functions) is needed. An elementary description holds however in the 
"logarithmic" case c(x) = 0. Then, (|1U|) and (|lip already tell that x_i and x are two roots of 
a(x) = 0. And as the polynomial a has degree 3 in TheoremCQ let a(x) = (x — x_i)(x — x )(x — (). 
Then, from (fl~3"D. 



-Cl- 



Ci 



xo 



■(4 



Xn—X 



{y-x-y[) ■ ■ ■ {y-i- y' n _ 1 )X 2 {x-i){x. 



xo 



(x_i- X n _ 2 ] 



(y-x- yx) ■ ■ ■ (y-x- y n -2)(x-i- x' ) • • • (x_ 



a(x' k ) (x fc - x_i)(x fc - x' ) 



c n y n (x) 



n-X 

2 a (x fe ) (x' k -x_i)(x' fc -x ) 
Ci 



-cr 



x 



x 



-(x n X n — 1) 



- y'l) ■ ■ ■ (y-x - y' n _ 1 )x 2 (x-i)(x- 1 - x ) • • • (x_ x - x„_ 2 ) 



yi 



n-X 



n 



(v- 



(x' k - C)(x 



Vk) 



y n - 2 )(x-i - x' ) ■ ■ ■ (x_i - x' n ) 



[ (xk - C)(x - y' k+1 )' 



(14) 
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3 A word on convergence 
3.1 Average behaviour 



We expect products occurring in (|13p or ()14[) to behave like powers, like 

n n 

n> - x k ) = no - + kk )) ~ $+(x) n - 



What is c &+(x) = expV+(x), where V+ is the complex potential of the distributions of x\P. 
For x' k , we write V-{x). For y, let us use the symbol W. 
The main behaviour of the n th term of (|14|) is therefore 



exp (n(W_(y_i) - + V+(x_i) - V_(s_i) + V_(C) 

-V + (C) + W+(x)-W_(x))). 



(15) 



Remark that we will only need V = V+ — V_ and W = W+ — W-. 

If is a general complex number, fill the whole complex plane and no convergence occurs. 

Let h be a real irrational multiple of a period u, then the same factors reappear approximately 
in the product after N steps if Nh is close to an integer times u. <l?(x) is the limit of the ./V th 
roots of such products. The various kh, for k = 1,2, ... ,N, modulo oj, fill uniformly the segment 
[0, u], and Xfc fill a curve which is the set of S(to + u),u E [0,u]: for any j in {1,2, . . . , N}, 
there is a A; such that kh is close to ju/N modulo oj. Indeed, let Nh be close to M^ui, with 
gcd(iV,Mjv) = 1. Then, 



kh 



M N 
N 



k + u 



kM N - j 
N ' 



to any j, there are integers k and m such that kM^ — mN = j (Bezout). 
So, we rearrange the product as 



$0) 



v 



Y[(x-S(ju/N + to)) 

3=1 



l/N 



exp 



log(x — + to)) du 



o 



As £ is the inversion of an elliptic integral of the first kind, 

" z dv 



u + 1 = 
we have 

3>(x) = exp 



log (a; — v) dv 



where {x n } is the locus = {£ (u + to)} 5 M £ [0, w]. The constant 1/lo is such that $>(a;) 
large x: 



x for 



dv 



So, let the complex potential 



V+(x) 



log (a; — v) dv 



'{xn} y/P{v) 
(V- will be used with x' n , and W± for y n and y' n ). 
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The formula for the potential will be linear after a convenient conformal map. 
One has the derivative 

w , , 1 f dv 
VAx) 



W J{x n } (X-V)y/P(v) 



with £ such that x = £ (£), dx/d^ = yJP{x). 

So, V+(x) and 1/_(x) are contour integrals on the loch filled by {x n } and {x' n } drawn by 
£(nh+to) and £(nh + t' n ). If x is between these two loch, the two contour integrals of dv , 

are the same for V' + (x) and V'_(x), up to the residue at v = x: 

., . . . , i / \ 2iri dV(x) 2iri 

ujy/P(x) of w 

We see that the real part of V remains constant on lines in the £-plane such that d^/uj is real, 
i.e., on parallel lines sharing the w-direction. 

Remember that the step h has been supposed to be a real multiple of co, so the arguments in 
arithmetic progression of step h in the £-plane of the elliptic functions defining a sequence x n , 
or y n , etc. happen to draw parallel lines with the w-direction! The real part of V(Q — V(x_i) 
occurring in fjl5[> is therefore 27r/|u;| times the distance between, say, if £ is the value of the 
elliptic function at and the line leading to the {x n } sequence. 

The remaining terms of Q15|) lead to a convergence behaviour dominated by 

exp(-n Im27r(£ a; - £c)/ w ); ( 16 ) 

where £ x is sent to x by the elliptic function. 

Of course, convergence holds while x is between the locus of x n and the corresponding locus 
(equipotential line) containing £. 

In a Jacobian setting, evaluation of (|16j) typically involves exp(—nirK'/K), well known in 
Zolotarev problems solutions and generalizations [8]. 

Rate of approximation has already been related to potential problems by Walsh [36^ chap- 
ter 9], in papers and books going back to the 1930s! See also Ganelius [?]• For more recent 
surveys and papers, the works of Bagby [3], and by Goncar and colleagues are recommended 

[a ei nni hue]. 

It is quite remarkable that configurations of particles in statistical physics [181 \19\ [20] are 
described in the same way than zeros and poles of rational approximations [3 [HI [91 [TU[ [TT] [T2"l 

[251 ES ES]. 



3.2 Exceptional cases 

The properties of the irrational number relating the step h to a period uj must also be conside- 
red [31]. Indeed, (fT4]) contains a division by a factor — y n -2) which is the difference of 
the values of a function of period to at arguments differing by an integer multiple of h, so that 
the result will be small whenever (n — l)h is close to an integer multiple of uj. The difference 
will never vanish, as h/ui is irrational, but could become VERY small infinitely often. The set 
of irrational h/u that could destroy the convergence estimate above is fortunately of vanishing 
measure in the set of real numbers, as shown by Hardy and Littlewood in [13] (and reproduced 
by Lubinsky in [2H pp. 854-855 and 871]). 
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